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The ground state spin- wave excitations and thermodynamic properties of two types of ferrimag- 
netic chains are investigated: the alternating spin-1/2 spin-5/2 chain and a similar chain with a 
spin-1/2 pendant attached to the spin-5/2 site. Results for magnetic susceptibility, magnetization 
and specific heat are obtained through the finite-temperature Lanczos method with the aim in 
describing available experimental data, as well as comparison with theoretical results from the semi- 
classical approximation and the low-temperature susceptibility expansion derived from Takahashi's 
modified spin-wave theory. In particular, we study in detail the temperature vs. magnetic field 
phase diagram of the spin-1/2 spin-5/2 chain, in which several low-temperature quantum phases are 
identified: the Luttinger Liquid phase, the ferrimagnetic plateau and the fully polarized one, and 
the respective quantum critical points and crossover lines. 

PACS numbers: 75.10.Pq, 75.10.Jm, 75.30.Kz, 75.40.Mg 
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I. INTRODUCTION 

Quasi-one-dimensional magnetic materials form a class 
of compounds with magnetic properties that above a 
characteristic temperature can be described through one- 
dimensional models [l[ . These include systems with a ro- 
tationally invariant singlet ground state (GS), modeled, 
for example, through spin-1 gapped and gapless critical 
spin-1/2 chains as wei l as more complex structures, 
such us ladders [l| and spin tubes Q . Typically, gap- 
less one-dimensional systems exhibit power-law decay of 
the correlation functions and can be understood through 
the Bethe ansatz 0] or the Luttinger Liquid theory [l|. 
In addition, in gapped one-dimensional (ID) systems the 
application of an external magnetic field B can suppress 
the gap and induce a quantum phase transition [J] to a 
Luttinger Liquid phase. In particular, an extensive study 
of the B - temperature (T) phase diagram of a spin-1/2 
gapped ladder system was recently carried out [H Q . 

Contrary to the above mentioned systems, quasi-lD 
ferrimagnetic compounds display GS spontaneous mag- 
netization and have ferromagnetic and antiferromagnetic 
(AF) spin-wave excitations. Usually the AF spin-wave 
mode is gapped and the magnetization curve exhibits 
a plateau, which can be explained by topological argu- 
ments 0. Ferrimagnetism can arise from the topology 
of the unit cell Q, as in the phosphates with chemical 
formula A 3 Cu3(P04)4, where A = Ca, Sr or Pb. These 
materials have three Cu 2+ spin-1/2 ions @ and can be 
modeled by a line of spin-1/2 trimer clusters QUEH with 
AF exchange couplings. Another class of ferrimagnets 
are mixed-spin compounds of type (A-X-B-X-)„, where 
A and B are two different magnetic components (single 
ions or more complex molecules) and X is a bridging lig- 
and. In particular, we are interested in compounds that 
can be modeled by spin-1/2 spin-5/2 chains (sS chains); 
this includes, for example, systems built from Mn 2+ and 
Cu 2+ ions linked through a dithioxalato ligand [HI, [l]| . 
Further, in the composition of some ferrimagnets, mag- 
netic elements can be organic radicals like the nitronyl ni- 



troxide free radicals (NITR), where R stands for an alkyl 
(methyl, ethyl) or aromatic group (phenyl). A family 
into this category consists of the Mn-NITR compounds 
[Til ] , for which there is an AF exchange coupling between 
the spin-1/2 radicals and the spin-5/2 Mn 2+ ions. 

In this work we present a numerical study of the 
thermodynamic properties of the ferrimagnetic chains 
illustrated in Fig. [TJ spin-1/2 spin-5/2 alternating 
chain (sS chain) and the spin-1/2 spin-5/2 alternat- 
ing chain with a spin-1/2 pendant attached to the 
spin-5/2 site (ssS chain). These chains can be re- 
spectively used to model the ferrimagnetic compounds 
CuMn(S 2 C 2 2 )2 • 7.5H 2 (denoted by CuMnDTO) [3 
and [Mn(NITIm)(NITImH)]C10 4 (denoted by MnNN) 
[l5l ]. whose crystal structures belong to the centrosym- 
metric monoclinic spac e g roup P2\jc [C^h)- The 3D or- 
dered phase observ ed Jl3l . ITa . Il6j ] in these compounds, 
and in similar ones [17] . at very low-T have been inten- 
sively investigated. In fact, magnetization measurements 
[TEL Il6j ] suggest that the canting of the ferrimagnetic mo- 
ments of the chains give rise to a 3D weak ferromag- 
netism below the critical temperature, although neutron- 
diffraction experiments [l7j in similar compounds in- 
dicate a canted AF structure. A common feature in 
these compounds is that the Mn 2+ ion has a 6 5 , 5/ 2 GS, 
thereby leading to a single ion anisotropy with no zero- 
field splitting in first order of perturbation theory. Very 
low-T magnetization measurements in a noncentrosym- 
metric orthorhombic compound [l7| . belonging to the 
space group P2i2i2i, suggest a single-ion anisotropy 
D/ks ~ 40 mK, which is much smaller than the intra- 
chain AF exchange couplings of the referred compounds 
[15. fl5l - [l7| . Therefore, a proper description of the 1D-3D 
magnetic transition may require, in general, anisotropic 
couplings, including the dipolar interaction. However, 
similarly to previous analysis [l~3|, [l5j], in order to de- 
scribe the ID ferrimagnetic properties of the compounds, 
we disregard anisotropy effects. 

The GS and the low-energy magnetic excitations are 
calculated through the Lanczos exact diagonalization 
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FIG. 1: (Color online) Schematic representation of the ex- 
change couplings and the GS long-range ferrimagnetic order- 
ing for the (a) sS and (b) ssS alternating chains. 



(ED) algorithm, while thermal properties are obtained 
by the finite-temperature Lanczos method (FTLM) [l8j]. 
We also explore the field-induced quantum phase transi- 
tions of these systems and discuss our results in light of 
experimental data, as well as predictions from the semi- 
classical approximation [l9j and the modified spin-wave 
(MSW) theory @. 

This work will unfold as follows: in Sec. II, we de- 
scribe the theoretical models and methods employed. In 
Sec. Ill, we estimate the model parameters suitable to 
describe the experimental data (susceptibility and mag- 
netization) of the related compounds. In Sec. IV, the one 
magnon bands and the specific heats of the two systems 
are presented and the main features discussed. In Sec. V 
we exhibit the T — B phase diagram of the sS chain and 
discuss in detail its quantum critical points and crossover 
lines, the Luttinger liquid phase and the plateau regions. 
In Sec. VI, we analyze the low-temperature behavior of 
the zero-field magnetic susceptibility and, finally, in Sec. 
VII we present a discussion of our relevant findings. 



The ssS chain with N c 
following Hamiltonian: 



[s,S,s'] cells is described by the 



H. 



ssS 



J2^s,-(S, + Sz + i) 

1=1 



1=1 



Sl - gf i B BS z , (2) 



where ,s = s' = 1/2 and S ~ 5/2, while J > and 
J' > 0. The Lieb and Mattis theorem assures that Sqs = 
N C \S - 2s| = 3V c /2, i. e., spin § per unit cell. The GS 
magnetic order of this chain is sketched in Fig. [lib). 

The FTLM [l8[ is based on the Lanczos diagonaliza- 
tion technique and random sampling. The fundamental 
relations used in FTLM for the calculation of an static 
quantity associated to an operator A are 
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where the sampling is carried over R random states |r), 
taken as initial states for a .M-step Lanczos procedure 
which results in M. approximate eigenvalues with re- 
spective eigenvectors \ifjj) in the 7V st -dimensional Hilbcrt 
space. The method allow us to calculate the tempera- 
ture dependence of the magnetization per unit cell m c , 
magnetic susceptibility per unit cell Xi and specific heat 

C through: m c = gfj, B ^-, x = g 2 ! mI <(S n^Jt } > and 

C = ^fc^rl^ i where ks is the Boltzmann constant. 
The total number of sites N = 2N C for the sS chain and 
N = 3N C for the ssS chain. In the computation we have 
used periodic boundary conditions, At = 50 for both 
chains and R = 40000 (50000) for the sS (ssS) chain. A 
full diagonalization study of the specific heat and suscep- 
tibility for the sS chain with N c = 3 can be found in Ref. 

m. 



II. MODELS AND METHODS 

The sS chain with a uniform exchange interaction 
J(> 0) and N c [s,S] cells is described by the following 
Hamiltonian: 



H sS = jJ2( S i + S *+i) ' s * - 9^bBS z 



(1) 



where s — 1/2 and S = 5/2, the (/-factor is assumed uni- 
form, hb is the Bohr magneton, B is an applied magnetic 
field in the z direction and S z is the operator for the z 
component of the total spin. This chain is bipartite, with 
N c sites with spin-5/2 in one sublattice and N c sites with 
spin- 1/2 in the other; therefore, the Lieb and Mattis the- 
orem [201 assures that the GS total spin, Sgs, is given 
by N C \S - s\ = 2N C , i. e., spin 2 per unit cell. The GS 
magnetic ordering of this chain is sketched in Fig. Ula). 



III. MAGNETIC SUSCEPTIBILITY AND 
MODEL PARAMETERS 



Through a semiclassical approach, in which the S spins 
are treated as classical variables, Seiden [lj| derived a 
closed formula for the magnetic susceptibility x- Li par- 
ticular, the quantity Tx{T) has a minimum at a temper- 
ature T m i n which is generally situated in a region where 
P J S < 1, a feature which has been known to be typi- 
cal of ID ferrimagnets. Similarly, a closed expression for 
the susceptibility of the ssS chain can also be established 

In Fig. [^a) we present data for the magnetic suscepti- 
bility of the compound CuMnDTO (from Refs. dlGS) 
together with FTLM, with J/k B = 44.8 K and g = 
1.90, and semiclassical-approximation (l9j results, with 
J/kB = 59.7 K and g = 1.9, for the sS chain. For 
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FIG. 2: (Color online) Molar susceptibility Xm times temper- 
ature T of the (a) sS and (b) ssS chains as a function of T. 
Experimental: data of (a) the compound CuMnDTO from 
Ref. Q]| and of (b) the MnNN compound from Ref. [H. 
FTLM: (a) N = 16, g = 1.90 and J = 4A.8K, (b) N = 18, 
g = 2.0, JJk B = 150 K and J'/k B = 255 K. Semiclassical: 

(a) Ref. [Til (with J/k B = 59.7 K, S = 2.5 and g = 1.9) and 

(b) Ref. QJ] (with J lk B = 141 K, .]' /k B = 250 K, S = 2.5 
and g = 2.0). 



the FTLM results, the estimation of J is made by us- 
ing the value of the minimum of the experimental curve: 



130 K. We sec that both the FTLM and the semi- 



classical approach agree with the experimental data in 
the mid- and high-temperature regimes. As the temper- 
ature is lowered below T m i n , \T increases and presents a 
maximum at T max = 7.5 K, which marks the onset of the 
tridimensional ordering. For a strictly one-dimensional 
system it is expected, from the Mermin- Wagner theo- 
rem [22| . that long-range order (LRO) may occur only at 
T = 0. Wc remark that for quantum ferromagnetic [23[ 
chains the correlation length diverges as 1/T and the sus- 
ceptibility as 1/T 2 ; further, the low-lying magnetic ex- 
citations of ferrimagnetic chains present a ferromagnetic 
character (see below) and the same referred critical be- 
havior is shown [24| to hold, which explains the increase 
in the curve of %T just below T = T m i n . 

In Fig. EJb) we present FTLM and semiclassical results 
fisl ] for the ssS chain, and experimental data of the 
MnNN compound. A profile similar to that of Fig. HJa) 
is observed with T m i n = 255 K. Taking g = 2.0, our 
estimative for the model parameters, J and J', are J' = 
1.7 J, with J/k B = 150 K and J'/k B = 255 K. Wc also 
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FIG. 3: (Color online) (a) Molar magnetization M m of the sS 
chain as a function of B at T — 4.2 K. Experimental data of 
the compound CuMnDTO from Ref. FTLM results for 

a chain with N = 16, J/k B = 44.8_K and g = 1.93. Brillouin 
paramagnet (BP) with total spin S = 15.0. (b) FTLM results 
for Mm of the ssS chain with N = 18 as a function of B at 
T = 15 K. BP with S = 8.5. 



note that the FTLM curve for this chain depart from the 
experimental one at a higher temperature than the curve 
for the sS chain. This behavior is in fact a finite size 
effect since the number of unit cells used in the FTLM 
calculation for the ssS chain (6 unit cells, with 18 sites) 
is effectively less than the number used for the sS chain 
(8 unit cells, with 16 sites). 

In Fig. |3]we compare our results at T = 4.2 K with ex- 
perimental data from Ref. (l3| . The temperature is lower 
than the one in which the maximum of the %T curve is 
observed, T max w 7.5 K. However, due to the low value 
of the coupling between chains ./inter-chain ~ 0.1 K, 
we expect that at T = 4.2 K and for fields higher than 
~ 0.1 T, the ferrimagnetic correlations along the chain 
are the relevant ones to determine the behavior of the 
magnetization as a function of B. Since the correlation 
length along the chains diverges 0, as 1/T, a finite 
number of unit cells are correlated at 4.2 K. Thus, we 
can treat the system as composed of independent lin- 
ear clusters, each cluster carrying a total spin S, and 
a superparamagnetic behavior is expected for the mag- 
netization curve. Within this context, we try to esti- 
mate the number of correlated cells in the chain from the 
experimental data shown in Fig. [SJb) by comparing it 
with the FTLM data and the molar magnetization of a 
Brillouin paramagnet (BP) with total spin S, given by 
M m (B, T) = Na9Hb{S- s)Bg[x), where N A is the Avo- 
gadro constant and B-g(x) is the Brillouin function. As 
shown in Fig. [21b), the experimental data is well de- 
scribed by the BP curve with S — 15, indicating that 
approximately 8 unit cells (size used in the FTLM calcu- 
lation) are ferrimagnetically correlated at 4.2 K. This en- 
forces the one-dimensional description of the experimen- 
tal magnetization for this temperature and field values, as 
well as the superparamagnetic behavior. Wc remark that 
the authors of Ref. [l3| estimate that approximately 10 
cells of the compound CuMnDTO are ferrimagnetically 
correlated at T — 7.9 K (just above T max ) for B = 0. 

For the ssS system, we find no published experimental 
data for the magnetization. However, considering the 
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FIG. 4: (Color online) ED results for the lower energy one- 
magnon bands, in units of magnetic field (using g = f .93), of 
the (a) sS and (b) ssS chains for the indicated values of N. 
Full lines are non- interacting spin- wave results from Ref. [2(|, 
while dashed lines are guide to the eyes. 



FTLM results, Fig. E^b), wc estimate that at T = 15 
K the number of fcrrimagnctically correlated unit cells is 
~ 6 (size used in the FTLM calculation) , due to the good 
agreement between the FTLM results and the BP curve 
with 5 = 8.75. 



IV. ONE-MAGNON BANDS AND SPECIFIC 
HEAT 

Due to the ferrimagnetic order of the GS, there are 
two kinds of elementary excitations in the systems: fer- 
romagnetic magnons, which lowers the total spin by one 
unit and AF magnons, which increases the total spin by 
one unit. The dispersion relations of the lower energy 
magnons are calculated, respectively, through 



u (q) 
u+(q) 



i{Sgs 
i(Sgs 



1,9) 
1,9) 



Egs 
Egs, 



(4) 
(5) 



where E m i„ (St , q) indicates the lowest energy in the 
total-spin sector St and lattice wavenumber q = 
2ttI/N c , with / = 0, 1,2, ...,N C - 1. 

In Figs. Hta) and@|b) wc display the lower energy one- 
magnon bands, in units of magnetic field, of the sS and 
ssS chains. For the sS chain, we also plot in Fig. QJa) 
non-interacting spin- wave (SW) results [2q |: 



J sw 



"SW 



(q) = -J(S-s)+uj q + gn B B, (6) 
(q) = J(S - s) +ui g - gfi B B, (7) 



where u q = J^J (S - s) 2 + 4s S sin 2 (q/2), s = 1/2 and 



5" = 5/2. We notice that in zero field the ferromag- 
netic excitation is gapless (which is expected from the 
spontaneously broken symmetry of the GS) and display 
a quadratic dispersion relation in the long wavelength 
limit, as predicted by conformal invariance [27j . while a 
gap A exists for the AF excitation. The ferromagnetic 
branch obtained through non-interacting SW theory for 
the sS chain is in good agreement with the ED data, 
while for the AF branch the value of the zero field gap 
A = 2 J(S — s) = 4 J departures from the ED value, as is 
often the case in other ferrimagnetic systems [28[ , due to 
quantum fluctuations effects. In fact, we estimate that 
in the thermodynamic limit A = 4.9046 J (3.88 J) for the 
sS (ssS) chain. 

In Fig. [SJa) we show the specific heat of the sS and ssS 
chains in zero field. Due to the LRO ferrimagnetic state 
at T = = B, with low-energy gapless ferromagnetic ex- 
citations, it is expected that C ~ VT. Another feature is 
the occurrence of double peaks [29Tj3l| ; it turns out that 
the main peak is well described by the Schottky formula 

= col^VW WhCTe 6 18 SchOUky 

gap and A is the amplitude parameter. The Schottky 
gap for the ssS chain 5 ss s ~ 4.1 J is in accord with the 
AF spin-wave gap A ss s ~ 3.9J. However, for the sS 
chain the value of 5 S $ (~ 3.4J) significantly departures 
from the AF spin- wave gap value: A s $ ~ 4.9 J, indicating 
strong influence of the lower-energy ferromagnetic exci- 
tations. Since these states have a total spin (St = S g — 1) 
lower than the one of the AF branch (St = S g + 1), we ex- 
pect that a field B can wash it out. In fact, the center of 
the AF branch [see Fig. QJa)] is found at A sS (0) = 5.4 J 
and is lowered in the presence of a magnetic field through 
A s s(B) = A s s(0) — gfisB. In Fig. Mjo) we present the 
specific heat for fields up to 103.8 T, and in Fig. (5Jc) we 
compare A S $(B) with the Schottky gap 5 s s(B). As we 
can see in the figure, the values of the two quantities are 
nearly equal for moderate values of B. 



V. T - B PHASE DIAGRAM 

The GS magnetization per unit cell, m c , of one- 
dimensional systems under an applied magnetic field can 
exhibit plateaus at values such that S c — m c = integer, 
where S c is the maximum total spin of a unit cell [7[. 
This condition implies that a plateau can be observed 
at values of m c differing from its saturated value by an 
integer number of spin flips. In particular, a magneti- 
zation plateau at 1/3 of the saturation magnetization 
was observed in the magnetization curve of the min- 
eral azurite [32|], which is gene rally modeled through 
the distorted diamond chain [331 ] . Other compounds ex- 
hibiting the 1/3 magnetization plateau are the trimer 
chain systems Cu2(P20eOH)2 [34[ and the phosphates 
A 3 Cu 3 (P0 4 )4, where A = Ca, Sr or Pb. Further, 
the thermal properties of a variety of models [35Tj37l | pre- 
senting plateaus in their magnetization curves were ana- 
lyzed in recent years and it was evidenced that the 1/3 
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FIG. 5: (Color online) FTLM results, (a) Specific heat of the 
sS and ssS chains for B — as a function of temperature T; 
full lines are the respective Schottky specific heats, (b) Molar 
specific heat C m of the sS chain for the indicated values of B 
(using g = 1.93); (c) B-dependent Schottky gap S s s(B) and 
spin- wave gap A s s(B). Dashed lines are guide to the eyes. 



magnetization plateau is also a characteristic feature of 
frustrated spin-S chains (38|. 

For the sS chain studied here, possible plateaus should 
be observed at m c = mpM = gf-i B (S — s) [Licb-Mattis 
(LM) magnetization] and m c = mpp = gpi B (S + s) 
(fully polarized magnetization), as confirmed by the nu- 
merical results shown in Fig. [6] (a). In zero field the 
GS is ferrimagnetic with gapless ferromagnetic excita- 
tions [Eq. ©], while for B ^ this mode acquires the 
gap Ap(B) = g/ipB. Also, as B increases, the gap 
for the AF mode [Eq. (JT])] decreases linearly with B: 
A AF (B) = A sS - g[i B B, and for g^pB > gnpB m = 
(A s s/2) ss 2. 45 J its gap is equal to the ferromagnetic 
one. At g^pB = A s s = gnE,B Ct AF the AF gap van- 
ishes and the system undergoes a quantum phase transi- 
tion (condensation of AF magnons, each carrying a spin 
+1) to a gapless Luttinger Liquid (LL) phase |4(J, with 
power-law decay of the transverse correlation functions. 
In fact, the quantum critical point B c ^af separates an 
incompressible phase (plateau) from a compressible one 
(LL phase). For B > B c ^af, & low-density of magnons 
is found in the system and the asymptotic singular form 
of the magnetization can be obtained fill ] by considering 
the system as a free Fermi gas or hard-core bosons. In 
this limit, the magnons will occupy single particle states 
with q — > and the dispersion relation, Eq. (JT)), can be 
used by replacing the linear spin- wave gap, Asw, by the 
computed gap in Fig. g|a), A sS = 4. 9046 J: 
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FIG. 6: (Color online) FTLM results for the sS chain, (a) 
GS magnetization per cell m c as a function of B for N = 16: 
full circles indicate the midpoints in the steps of the magneti- 
zation of the finite-size system and edges of thermodynamic- 
limit plateaus [39| . while colored full lines indicate the results 
from the free-fermion model (see text), (b) m c (curves with 
symbols) and \ (full curves) as a function of B for the listed 
values of T; B m is indicated by the red diamond. 



g[i B [B — B c ^af) and Asw = 2J(5 — s). The energy 
density can thus be written (fermionic map) as 



dk 
2^ 



(9) 



where = v k /2A s s, kp = 7m, and n is the density 
of particles. The value of n for a prescribed fi can be 
obtained from the condition cLe = 0: 



- 2B c ,af r- 



(10) 



which implies 



= 2 



SUB 



J 



g^B ^B C ,AF 



5tt 2 



S/B - B C ,AF. (11) 



In Fig. EJa) we show the very good agreement between 
the numerical data and m c given by Eq. (|11[) in the LL 
Phase. 
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FIG. 7: (Color online) FTLM results for the low-T phase 
diagram of the sS chain: contour plot indicates the magneti- 
zation per cell m c . The critical point at B — (black dia- 
mond), the inflection point of m c at g/igB = giisBm = A/2 
(red diamond), the quantum-critical points (black triangles) 
at B = B c ,af and B — B c ,fp, crossover lines (white circles) 
and their asymptotic behavior (full lines) are also indicated. 



The gapless LL phase ends at the quantum critical 
point B = B c pp: the system becomes fully polarized 
(FP) and presents gapped low energy excitations. The 
two one-magnon excitations from the FP state, both car- 
rying a spin -1, can be exactly obtained [42| and the lower 
one has a dispersion relation given by 

lu fp = -J(s+S)- Jy/(S- s) 2 +4sScos 2 (q/2)+gfji B B, 

(12) 

which implies g^pBc^FP = 2J(s + S) = 6J, in accord 
with the numerical results [Fig. EJa)]. For B < B c pp, a 
low density of magnons is observed in the system and 
the same arguments used to obtain Eq. (fTOj) can be 
used in this case. For q — > 0, Eq. (|T2"j) can bc writ- 
ten as Eq. (HJ with v = jV2sS = Jy/b/2, A F = 
2J(s + S) = 6 J and fi = gfiB{B c ,FP — B), which implies, 
from Eq. (|10j) . that the density of magnons is given by 

n — 9j±s. \J iB ^ p \j B c ,pp — B, and m c now reads: 

= 3 rV 5 2 V B c,fp-B, (13 

which is plotted in Fig. H2a) and is also in very good 
agreement with the numerical data. 

In Fig. EJb) we present FTLM data for m c and x vs. i? 
for T/0. We first notice that the magnetization in zero 
field is null and the system is in the thermal paramagnetic 
state, as expected from the Mermin- Wagner theorem 
[22j . Increasing B in the low-temperature regime, the LM 
(or ferrimagnetic) plateau is exponentially reached and 
the magnetization exhibits an inflection point at B = B m 
[red diamond in Fig. H>(b)] that marks the changing of 



the gapped low-energy excitations from ferromagnetic 
(B < B m ) to AF magnons (B > B m ): the ferromag- 
netic (antiferromagnetic) magnons are exponentially ac- 
tivated and m c is lower (higher) than ttilm = 9(J-b(S—s). 
Also, by the same token, the FP plateau is exponentially 
reached from below for fields higher than B c pp. Fur- 
thermore, the singular form of the magnetization near 
the quantum critical points (B = B c af and B = B Ct pp 1 
with T = 0), which implies x — > oo, are thermally 
smoothed out and the singularities in the susceptibility 
evolve into local maxima, thus providing the determina- 
tion of the crossover lines. The LL phase, with linear 
dispersion relation ~ q } is expected [40( between the two 
local maxima for a given T [see, e. g., the susceptibility 
curves for k B T = 0.10J and 0.20J in Fig. |B(b)] with 
the two local maxima indicating a crossover to a region 
in which the excitations follow a non-relativistic disper- 
sion relation ~ q 2 , as previously discussed. On the other 
hand, as T increases, the LL phase ends and a single 
maximum is observed in the susceptibility curves (see, 
c. g., the susceptibility for ksT = 0.40J). This single 
maximum defines a crossover from the regime in which 
the physics is determined by the excitations from the LM 
plateau to a regime in which the FP plateau is the rele- 
vant one. For sufficiently high temperatures, the system 
looses all information about the T = LM magnetiza- 
tion plateau and the effect of B is to bring the system 
from the thermal paramagnetic state to the FP state at 
higher magnetic fields (see the case fegT = 2.00J). 

In Fig. [7] we present the contour plot of m c in the 
T — B plane and a schematic phase diagram. The T — B 
crossover lines enclosing the region of the LL phase, lim- 
ited at T = by B = B c ^af and B = B Ct pp, are obtained 
[43( 1 from the local extrema of m c (T) vs. T for a given 
B, as shown in Fig. Ufa). Further, as B — > B c these 
crossover lines follow a universal function 01 : a\B — B c \ 
with a = 0.76238; as shown in Fig. our numerical data 
confirm this asymptotic behavior for the two quantum 
critical points at B = B c ,af and B = B c F p. More- 
over, as T increases beyond the crossover lines of the 
two plateaus, gapless phases are reached [4(| H |. In 
addition, by increasing B under a fixed T, local max- 
ima are observed in the specific heat C(B) per spin, as 
displayed in Fig. [Sfb). These features are used to es- 
timate Q the crossover lines related to the LM plateau 
and FP plateau shown in Fig. [7j in particular, we notice 
that T ~ \B — B c \ as the lines reach the corresponding 
quantum critical points 0, E3|. Last, we stress that the 
crossover lines and the LL instability lines meet at the 
quantum critical points, thus delimiting the respective 
quantum critical region [4(| l43l - |45j ; in each region the 
system is thus governed by the quantum critical point 
with dynamical exponent z — 2 associated with the ex- 
cited magnons, as discussed above. On the other hand, 
the magnon densities [40l |. n, given by (m c /g/is) — 2 
and 3 — (m c l g[ip) for the quantum critical point at 
B = B c ^af and B = B c ,pp, respectively, follow a univer- 
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FIG. 8: (Color online) FTLM results for the sS chain, (a) 
Magnetization per cell m c as a function of T for fields qhbB/ J 
from 4.85 to 5.95 in steps of 0.05 (from below to top); trian- 
gles indicate local maxima associated with the LL crossover 
lines (see text), (b) Specific heat per cell C as a function 
of B for the indicated values of temperature. Scaling of the 
magnon density n around the quantum critical points at (c) 
B = B c ,af and (d) B = B C , FP . 



sal function of T and \B — B c \/T: 



k B T 
J 



f 



\B-B C 
T 



(14) 



as shown in Figs. |5] (c) and [5] (d). A better scaling 
behavior is observed for B < B c ,af {B > B c pp) in Fig. 
M (c) [H (d)] since for B > B C . AF (B < B C . FP ) the zone 
of influence of the quantum critical point at B = B c pp 
(B = B c _af) merges with the zone of influence of the 
point at B = B c ,af (B = B Ci fp)- The guideline ksT = 
giiBB in Fig. [7] is discussed below. 

In Fig. [9] we present the contour plot of C/T in the 
T — B phase diagram (BJ, including the above-discussed 
crossover lines. At the plateaus, C/T ^ as T 
due to the gaps, as evidenced in the plot. As we can 




CN . 
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FIG. 9: (Color online) FTLM results for the low-T phase 
diagram of the sS chain: contour plot indicates C/T. The 
critical points and crossover lines are indicated as in Fig. [7] 



see, the guideline ksT = gusB do not coincide with 
the local maxima of C(B) in the low-B region [see Fig. 
[5{b)] due to the LRO ferrimagnetic state at T = = 
B: since C - y/T, C/T ->ooasT^0atB = 
and an enhancement in the intensity of C/T is observed 
near T = = B. In spite of this fact, the plot shows 
a depression in the values of C/T near the T = LM 
plateau which, by increasing T, varies in a symmetrical 
fashion with respect to B = B m (dome-shaped) and is 
limited by the ksT = g^BB and fcsT = g\iB\B c ,AF — B\ 
asymptotic crossover lines. Further, the LL dome is also 
clearly seen and the crossover lines of the FP and LM 
gapped phases can be visualized. 

Next, we exhibit in Fig. [TU]the magnetization of the 
ssS chain at T — 0. For this chain, the first plateau is 
found at m c = g^B(S — 2s), i. e., the LM plateau, and 
the second is the FP plateau at m c = g/is (S+2s); the LL 
phase is expected to occur between these two plateaus. A 
third plateau could be found at m c = g^BS; however 
our numerical shows no evidence of this plateau. We 
remark that we did not perform a detailed analysis of 
the T — B phase diagram of this chain, but we expect 
that it should display similar features already reported 
for the sS chain. 

The huge values of the quantum critical magnetic fields 
of the CuMnDTO (sS chain) and MnNN (ssS chain) com- 
pounds, make the experimental investigation of the full 
T—B phase diagram of these systems very difficult. How- 
ever, magnetic phase transitions induced by very large 
magnetic fields (up to 400 T) in the low-temperature 
regime have been reported [4o| . Further, materials phys- 
ically described by similar models may have lower values 
for the exchange coupling and thus a more experimen- 
tally accessible phase diagram. 

We also mention that fcrrimagnctism can be destabi- 
lized by competing (or frustrating) interactions [U |48| , 
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FIG. 10: (Color online) FTLM results for the magnetization 
per cell m c of the ssS chain as a function of field B at T = 
for TV = 18. Full circles indicate the midpoints in the steps 
of the magnetization of the finite-size system and edges of 
thermodynamic- limit plateaus [3^ |. 
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FIG. 11: (Color online) X T 2 vs. T for the spin-1/2 ferromag- 
netic chain in the low and very low (IcbT « J) temperature 
regions. FTLM results for N = 8 and TV = 24. Fitting of the 
FTLM results for N = 8 using [i + a (^)2 + ai(^)] 
(see text). MSW results up to second order in ksT / J . The 
inset shows the results for ksT < 0.1J. 



which can give rise to other critical points. Unconven- 
tional ferrimagnetism (non-bipartite lattices) was indeed 
found in one-dimensional frustrated structures [4QL l50j 
and in the Kagome lattice [El]]. Further, the mag- 
netocaloric effect in the kinetically frustrated diamond 
chain was recently investigated |52j ■ 



VI. LOW-TEMPERATURE MAGNETIC 
SUSCEPTIBILITY 

We now consider the temperature regime where ferro- 
magnetic excitations tend to be a predominant feature. 
In order to test and illustrate the accuracy of the FTLM 
in describing the susceptibility behavior at very low tem- 
peratures, we have calculated the susceptibility of the 
spin-1/2 linear ferromagnetic chain; the results for xT" 1 
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FIG. 12: (Color online) \T 2 vs. T for the sS chain. Exper- 
imental data of the compound CuMnDTO from Ref. [la ]. 
FTLM results for TV = 16. Fitting of the FTLM results for 
TV = 16 using [§ + a (^)i + oi(^)] (see textL MSW 
results up to second order in UbT I J from Ref. [2Sj]. The 
semiclassical results are from Ref. [la ]. 



as a function of T are shown in Fig. [TT] for systems with 
8 and 24 sites. The crossover to zero of the FTLM results 
as T — > is due to finite size effects. We note that for 
(ksT/J) > 0.3 the curves for the two chain sizes super- 
impose, thus suggesting that the thermodynamic-limit 
behavior has already been within numerical accuracy. 
Also, in the temperature range 0.06 < (fc^T/J) w 0.1, 
the results for the larger system is in good agreement 
with the expansion formula from Takahashi's MSW the- 
ory [23|], which up to second order in t = ksT/J reads: 
_xJ_ _ t -2 



obtain 

xJ 



3* 



-22 S %At2 + sA 2 t + O [ti 
0.582597 and g = 2. For s 



, where 
1/2 we 



(g^B ) 2 



— + 0.145649P + 0.16971* + 0(t 



(15) 

We stress that in the range < (k B T/J) < 0.1, Eq. 
(|15l) is in very good agreement with predictions from 
the Bethe-ansatz approach [53|, while the fitting of the 
FTLM results for N = 8 and 0.5 < (k B T/J) < 0.9, yields 
do = 0.140 and a\ — 0.186, in good agreement with the 
MSW coefficients. 

We now turn our attention to the low-temperature 
regime of the sS-chain susceptibility displayed in Fig. [T2] 
Firstly, we note that for (k B T/J) > 0.5 the FTLM results 
for TV — 14 (not shown) and N = 16 (8 cells) coincide, in- 
dicating that the thermodynamic limit has been attained 
in this temperature range. The experimental data nor- 
malized by J/k B = 44.8 K (g = 1.88) and J/k B = 59.7 K 
(g = 1.9) show the expected agreement with the FTLM 
results and the semiclassical formula, respectively, as al- 
ready displayed in Fig. dta). The MSW results comes 
from the expansion formula derived by Yamamoto et al. 
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|29|, which up to second order in t reads: 



Ss(S - s) 2 



(Ss)2(S-s)iAt2 



+(S - s)AH + O (* = 



(16) 



A relevant aspect of this expansion is that for S = 2 s 
we recover the Takahashi expansion for the ferromagnetic 
linear chain of spin s, which reinforces that the ferromag- 
netic excitation is the relevant one at low temperatures 
[H[. Setting s = 1/2 and S = 5/2 in Eq. US we obtain 



xJ 



= t~ 



- + 1.842334*2 + 0.678840* + O it * 



(17) 

The FTLM results can be fitted by a function of the 
form [§ + a (^5^)5 + ai(^I)]. Guided by our studies 
on the spin- 1/2 ferromagnetic chain, we have chosen the 
interval 0.5 < (fcsT/J) < 0.9 to fix the values of do and 
ail &o = 1-28 and a\ = 0.69, which can be compared 
with those in Eq. (fT7|) , and implies a good agreement for 
the integer-power coefficient and an order-of-magnitude 
agreement for the half-integer power coefficient. 

One should notice that the FTLM results and the ex- 
perimental data crossover to zero as (ksT/J) —> 0, in- 
stead of approaching the constant value sS(S — s) 2 /3 = 
5/3. Here one must distinguish two effects: with re- 
spect to FTLM, this is evidently a manifestation of finite- 
size effects, while for the experimental data one can at- 
tribute this to the 1D/3D crossover that takes place be- 
low T — 7.5 K (see Sec. III). In fact, in the 3D region 
the susceptibility behaves as \ ~ T -7 , with the critical 
exponent 7 < 2, implying that %T 2 — > as T — » 0. 



dimensional models related to quasi-one-dimcnsional fer- 
rimagnetic compounds: CuMnDTO and MnNN. In fact, 
the models are associated to two types of ferrimagnctic 
chains: the alternating spin- 1/2 spin-5/2 chain and the 
spin-1/2 spin-5/2 alternating chain with a spin-1/2 pen- 
dant attached to the spin-5/2 site. The finite temper- 
ature Lanczos method proved quite reliable, except at 
very low temperatures where finite-size effects hinder its 
accuracy. A particular feature of these systems is the 
presence of gapless ferromagnetic and gapped AF spin- 
wave (magnon) branches in zero field. As the magnetic 
field is increased, the low-energy excitation changes from 
ferromagnetic to AF and the magnetic field vs. tem- 
perature phase diagram displays characteristic crossover 
lines which distinguish these systems from spin-1 Hal- 
dane chains and two-leg ladder models. In particular, 
for the sS chain we have identified the quantum criti- 
cal points and the crossover lines, the Luttinger liquid 
phase, the ferrimagnetic (LM) and the fully polarized 
plateaus. The values of the exchange coupling param- 
eters of the compounds discussed in the text are in- 
deed very high. However, magnetic phase transitions 
induced by very large magnetic fields (up to 400 T) in 
the low-temperature regime have been experimentally 
investigated [4q[. Also, other compounds described by 
similar models can have lower values of the exchange 
parameters and a more experimentally accessible phase 
diagram. We expect that this work stimulates experi- 
mental and theoretical research with focus on the phase 
transitions induced by an applied magnetic field in the 
low-temperature regime of the large class of quasi-one- 
dimensional ferrimagnetic compounds. 
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